Quantum mechanical evolution towards thermal equilibrium 



00 

o 
o 

(N 
o 

Q 

(N 



^ ■ 

^3 : 



> 

in 

00 

m 
oo 

o 



X 



Noah Linden/ Sandu Popescu/ '' Anthony J. Short/ and Andreas Winter^' ^ 

^Department of Mathematics, University of Bristol, University Walk, Bristol BS8 ITW, U.K. 
^H.H.Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 ITL, U.K. 
^Hewlett-Packard Laboratories, Stoke Gifford, Bristol BS12 6QZ, U.K. 
^Department of Applied Mathematics and Theoretical Physics, University of Cambridge, 
Centre for Mathematical Sciences, Wilberforce Road, Cambridge CB3 OWA, U.K. 
^Centre for Quantum Technologies, National University of Singapore, 2 Science Drive 3, Singapore 117542 

(Dated: 12 December 2008) 

The circumstances under wtiich a system reaches thermal equilibrium, and how to derive this from 
basic dynamical laws, has been a major question from the very beginning of thermodynamics and 
statistical mechanics. Despite considerable progress, it remains an open problem. Motivated by this 
issue, we address the more general question of equilibration. We prove, with virtually full general- 
ity, that reaching equilibrium is a universal property of quantum systems: Almost any subsystem in 
interaction with a large enough bath will reach an equilibrium state and remain close to it for almost 
all times. We also prove several general results about other aspects of thermalisation besides equi- 
libration, for example, that the equilibrium state does not depend on the detailed micro-state of the 
bath. 



Leave your hot cup of coffee or cold beer alone for a 
while and they soon lose their appeal - the coffee cools 
down and the beer warms up and they both reach room 
temperature. And it is not only coffee and beer - reach- 
ing thermal equilibrium is a ubiquitous phenomenon: 
everything does it. Thermalization is one of the most 
fundamental facts of nature. 

But how exactly does thermalization occur? How can 
one derive the existence of this phenomenon from the 
basic dynamical laws of nature (such as Newton's or 
Schrodinger 's equations)? These have been open ques- 
tions since the very beginning of statistical mechanics 
more than a century and a half ago. 

One - but by no means the only - stumbling block 
has been the fact that the basic postulates of statistical 
mechanics rely on subjective lack of knowledge and en- 
semble averages, which is very controversial as a phys- 
ical principle. Recently however there has been signif- 
icant progress: it was realized that ensemble averages 
and subjective ignorance are not needed, because indi- 
vidual quantum states of systems can exhibit statistical 
properties. This is a purely quantum phenomenon, and 
the key is entanglement, which leads to objective lack of 
knowledge. Namely, in quantum mechanics even when 
we have complete knowledge of the state of a system, 
i.e. it is in a pure state and has zero entropy, the state of 
a subsystem may be mixed and have non-zero entropy. 
In this situation we cannot describe the subsystem by 
any particular pure state, and for all purposes it behaves 
as if we have a lack of knowledge about what its pure 
state is (i.e. it behaves as a probability distribution over 
pure states). This is in stark contrast to classical physics 
where complete knowledge of the state of the whole sys- 
tem implies complete knowledge of the state of any sub- 
system, and hence probabilities can only arise as purely 
subjective lack of knowledge (i.e. the subsystem has a 
well-defined state only we don't know what that is). 

This approach has become a very fruitful direction of 



research in recent years Ht S |3|], see also irnportant ear- 
lier work L4,j5,^iZil and numerical studies la]. 

Most notably, it was shown in (it H, |3 that almost all 
(pure) states of a large system are such that any small 
subsystem is in a canonical state, as if it were the result 
of averaging (with equal probability) over all possible 
states respecting the desired macroscopic constraints. 
However, the above result is limited in that it refers to 
the state at a given time, and deals only with the case of 
generic states. But what we are most interested in is not 
the situation at one given moment of time but in the time 
evolution - in particular, under what circumstances sys- 
tems reach equilibrium, and how much they fluctuate. 
Furthermore, the states of most interest, namely those 
initially far from equilibrium, are non-generic. These 
dynamical aspects are the subject of the present paper. 

Thermalisation seems a very straightforward process 
- put a system in contact with a large enough thermal 
bath and the system will reach equilibrium at the same 
temperature. Closer analysis reveals however that the 
process of thermalisation actually contains many differ- 
ent aspects and we can decompose it into the following 
elements. They constitute a rough roadmap of what one 
has to show in order to demonstrate that a particular 
physical system thermalises. 

1. Equilibration. We say that a system equilibrates 
if its state evolves towards some particular state 
(in general mixed) and remains in that state (or 
close to it) for almost all times. As far as equili- 
bration is concerned, it is irrelevant what the equi- 
librium state actually is. In particular, it need not 
be a thermal state (such as a Boltzmannian distri- 
bution) and indeed may depend on the initial state 
of the subsystem and / or on the initial state of the 
rest of the system (i.e. bath) in an arbitrary way. 

2. Bath state independence. The equilibrium state of 
the system should not depend on the precise ini- 
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tial state of the bath. That is, when defining a bath, 
we should only specify its macroscopic parame- 
ters, e.g. its temperature; when the system reaches 
equilibrium, the equilibrium state should depend 
only on the temperature of the bath. 

3. Subsystem state independence. If the subsystem is 
small compared to the bath, the equilibrium state 
of the subsystem should be independent of its ini- 
tial state. 

4. Boltzmann form of the equilibrium state. Under cer- 
tain additional conditions on the Hamiltonian (es- 
pecially the interaction term) and on the initial 
state, the equilibrium state of the subsystem can 
be written in the familiar Boltzmannian form ps — 

Realizing that thermalization can be decomposed in 
this way has major consequences. Firstly, it allows us to 
address each aspect separately. Secondly, and more im- 
portant, it allows us to greatly expand the scope of our 
study. Indeed, we will consider equilibration as a gen- 
eral quantum phenomenon that may occur in situations 
other than those usually associated with thermalization. 
In particular we need not restrict ourselves to: standard 
thermal baths (that are described by a given tempera- 
ture or restricted energy range), weak or short range in- 
teractions between the system and the bath, Boltzman- 
nian distributions, situations in which energy is an ex- 
tensive quantity, etc. Furthermore, we can consider sit- 
uations in which the systems do not reach equilibrium, 
and prove results about the bath or subsystem indepen- 
dence properties of the time-averaged state. 

In this paper, with very weak assumptions, we prove 
the first two elements above - equilibration and bath 
state independence. That is, with virtually full gener- 
ality, we prove that reaching equilibrium is a universal 
property of quantum systems, and that the equilibrium 
state does not depend on the precise details of the bath 
state. After introducing the setup and basic definitions 
in the next section, we go through the above programme 
step-by-step. 



1. SETUP AND DEFINITIONS 

In this section, we describe our general setup, and in- 
troduce the basic notations and definitions that will be 
used throughout the paper. 

The system. We consider a large quantum system, de- 
scribed by a Hilbert space H. We decompose this system 
into two parts, a small subsystem S, and the rest of the 
system that we refer to as the bath B. Correspondingly, 
we decompose the total Hilbert space asH = Hs ^Hb, 
where Hs and Hb (of dimensions ds and ds) are the 
Hilbert spaces of the subsystem and bath respectively. If 



either part is infinite-dimensional, we introduce a high- 
energy cut-off to render its dimension finite (and elimi- 
nate interaction terms from the Hamiltonian that would 
take the state outside the allowed subspace). Note, how- 
ever, that at this stage we do not ascribe the subsystem 
or bath any special properties. The subsystem S could be 
a single particle, a cluster of particles next to each other, 
or a number of isolated particles spread throughout the 
bath. It may even be something far more abstract, not 
necessarily related to particles at all: Any arbitrary de- 
composition of the Hilbert space into a tensor product 
{H ~ 'Hs®'Hb) actually defines a subsystem and a bath. 

The Hamiltonian.The evolution of the total system is 
governed by a Hamiltonian 

H = Y,Ek\Eu){Eul (1) 

k 

where \Ek) is the energy eigenstate with energy Ek- 
Throughout this paper, we consider the Hamiltonian 
to be completely general, except for the following con- 
straint: That it has non-degenerate energy gaps. 

We say that a Hamiltonian has non-degenerate en- 
ergy gaps if any nonzero difference of eigenenergies 
determines the two energy values involved. I.e., for 
any four eigenstates with energies Ek, Ec, Em and En, 
Ek ~ Ee = Em — En implies k = £ and m — n, or k — m 
and £ = n. Note that this also implies that the energy 
levels are non-degenerate. 

An important physical implication of this assump- 
tion is that the Hamiltonian is fully interactive, in the 
sense that no matter how we partition the system into 
a subsystem and bath they interact. Indeed any non- 
interacting Hamiltonian H = Hg + Hb has multiple de- 
generate energy gaps. Note in particular that for non- 
interacting Hamiltonians the energy is a sum of the sys- 
tem energy and bath energy. Any four energies satisfy- 
ing El = Ef + Ef, E2 = Ef + E^, E3^Ef+ E^ , and 
i?4 = E2 + E2 will lead to a degenerate gap. 

We emphasize that the restriction to Hamiltonians 
that have no degenerate energy gaps is an extremely 
natural and weak restriction. Indeed, adding an arbi- 
trarily small random perturbation to any Hamiltonian 
will remove all degeneracies (although such changes 
may take a long time to significantly influence the evo- 
lution of a state, here we are not concerned with the time 
scales). 

Note that except for the Hamiltonian having non- 
degenerate energy gaps, we allow it to be completely 
general. In particular, in a system composed of particles, 
it need not describe nearest-neighbour or bipartite inter- 
actions, but could contain interactions between all parti- 
cles simultaneously. Hence energy need not be an (even 
approximately) extensive quantity, as normally consid- 
ered in statistical mechanics. 

Notation. We denote by \^{t)) the global pure state of 
the system (i.e. of the subsystem and bath together) at 
time t. It is also convenient to write this state as a density 
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matrix p{t) = The state of the subsystem 

Ps{t) is obtained by tracing the global state p{t) over the 
bath. I.e. ps{t) — Tr^ p(t). Similarly the state of the bath 
at time t is described by psit) = Trs p{t). 

It is useful to define the time-averaged state of the sys- 
tem oj, which is given by 

CO = {p{t))^ = lim - f p(t)dt. (2) 

Similarly we define ujs and ojb as the time-averaged 
states of the system and bath respectively. 

It is also convenient to introduce the notion of the ef- 
fective dimension of a (mixed) state p by 

This tells us, in a certain sense, how many pure states 
contribute appreciably to the mixture. In particular a 
mixture of n orthogonal states with equal probability 
has effective dimension n. Unlike the support of the 
density matrix, this notion captures the probabilistic 
weight of different states in the mixture, and is continu- 
ous. 

It will also be important in what follows to consider 
the distance between two density matrices. There are 
many possible ways to define such a distance; here 
we use a strong and very natural distance, the trace- 
distance D{pi,p2) = \ Tr y/ {pi — P2Y- The trace- 
distance characterises how hard it is to distinguish 
two states experimentally (even given perfect measure- 
ments). When it is small, the two states are effectively 
indistinguishable. More precisely, it is equal to the max- 
imum difference in probability for any outcome of any 
measurement performed on the two states. Further- 
more, the maximum difference in the expectation val- 
ues of any operator A on the two states is the range of 
the eigenvalues of the operator times the trace-distance 

Finally, by Pr*{ } we denote the proportion of states 
in a Hilbert space having a particular property, accord- 
ing to the natural Haar measure. 

2. EQUILIBRATION 

We now come to the central result of our paper: Ev- 
ery pure state of a large quantum system that is composed 
of a large number of energy eigenstates and which evolves 
under any arbitrary Hamiltonian (with non-degenerate en- 
ergy gaps) is such that every small subsystem will equili- 
brate. That is, every small subsystem will evolve to- 
wards some particular state (in general mixed) and re- 
main in that state (or close to it) for almost all times. 

To understand the requirement that the state contains 
a large number of energy eigenstates we first note that 
this means that there will be a lot of change during 
the time evolution. Indeed, a single energy eigenstate 



does not change at all. (In this trivial sense, every en- 
ergy eigenstate is already equilibrated. We however 
are looking for systems that do not start in equilibrium 
but evolve towards it.) But why do we need signifi- 
cant change during the time evolution? In order for 
equilibration to occur, we need that some information 
about the initial state of the subsystem leaves the sub- 
system. Indeed, suppose we start with a subsystem far 
from equilibrium. Then it will evolve through a num- 
ber of distinct states on its path towards equilibrium. 
This implies that the state of the whole system (i.e. sub- 
system and bath) will also evolve through a number of 
distinct states. However, suppose now that the subsys- 
tem reaches equilibrium. By definition, this means that 
its state doesn't change (significantly) any more. Nev- 
ertheless, due to unitarity the state of the whole system 
cannot simply stop evolving, and it must continue to go 
through distinct states at the same rate as initially. In 
order for recurrences of the non-equilibrium state of the 
subsystem to occur very infrequently, the states of the 
whole system in which the subsystem is far from equi- 
librium must be only a small fraction of the total num- 
ber of distinct states through which the whole system 
evolves. Hence the whole system must evolve through 
very many states. In order for this to be possible, the 
state of the total system must be composed of many en- 
ergy eigenstates. 

We have argued above that for a subsystem initially 
out of equilibrium to reach equilibrium, the state of the 
whole system must go through many distinct states. 
What we prove in the present paper is that this is also 
sufficient. That is: 

Whenever the state of the whole system goes 
through many distinct states any small subsys- 
tem reaches equilibrium. 

Another way of looking at the time evolution of the 
state is to investigate what happens to the state of the 
bath. As mentioned above, due to unitarity, the state of 
the whole system must continue evolving at the same 
rate even when the state of the subsystem is at equi- 
librium and doesn't change. The total state can then 
change in two ways: through changes in the correla- 
tions between the subsystem and the bath and through 
changes of the state of the bath itself. What we prove is: 

Whenever the state of the bath goes through 
many distinct states, any small subsystem 
reaches equilibrium. 

Furthermore, we use our proof of this (and the earlier 
indented statement) to show that equilibration occurs in 
initial product states of the subsystem and bath, for al- 
most all initial states of the bath. 

The notion of evolving through many different dis- 
tinct states is mathematically encapsulated by the effec- 
tive dimension of the time-averaged state, d {uj) where 
OJ = {p{t))^. The connection between this and the num- 
ber of energy eigenstates is easily seen by expanding 
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\^{t)} (setting h= 1 for convenience) as 

|vl/(i))=^c,e-*^'^-*|i?fc), (4) 

k 

where X^fc kfcP = 1/ and hence 

k,l 

From the non-degeneracy of the energy levels (implied 
by non-degenerate energy gaps), 

u; = {p{t))^^J2\^,\'\^'^){E,\, (6) 

k 

giving 



Similarly, the notion of evolving through many differ- 
ent distinct states of the bath is mathematically encap- 
sulated by the effective dimension of the time-averaged 
state of the bath, d'^^{u!B) where ujb = {pB{t))f Note 
that in general, we expect the time-evolved state to ex- 
plore a significant portion of the bath state space, which 
is much larger than the subsystem's state space. Hence 
we anticipate that d °* (tJs ) will be much larger than ds ■ 

We are now in the position to formulate our first theo- 
rem in rigorous mathematical terms. A quantity central 
to this result is D{ps{t), ujs), the distance between ps{t), 
the state of the subsystem at time t, and its time aver- 
age, LOs = {ps{t))f In general ps{t) fluctuates around its 
time average ujs and the distance between them changes 
over time. To characterize these fluctuations we look at 
the time average of the distance {D{ps{t),uJs))f When 
this average is small, the subsystem must spend almost 
all of its time very close to ujs- In other words, when 
{D{ps{t),u}s))t is small, the system equilibrates to ojs- 

Theorem 1 Consider any state \^{t)) e Ti. evolving under 
a Hamiltonian with non-degenerate energy gaps. Then the 
average distance between ps{t) and its time average ujs is 
bounded by 

The proof of this result is given in appendix lAl 

By bounding {D{ps{t) , i^s)) f theorem tells us that 
the subsystem will equilibrate whenever the effective di- 
mension explored by the bath d°^{ujB) is much larger 
than the subsystem dimension ds, or whenever the ef- 
fective dimension explored by the total state d'^^{uj) is 
much larger than two copies of the subsystem (dimen- 
sion (i|). In other words, if {D{ps{t), 0Js))t is small, the 
system spends almost all its time close to the equilib- 
rium state. Indeed, as distances are always positive, it 



is easy to see that the proportion of the time for which 
D{ps{t), uJs) is more than K times {D{ps{t), t^s))j (with 
K an arbitrary positive constant) must be less than 1/K. 
If we assume that the energy eigenvalues of H have no 
rational dependencies (which is much stronger than the 
non-degenerate energy gaps condition, but still one that 
holds for generic perturbations of the Hamiltonian) we 
can improve the bound on the proportion of time the 
subsystem spends away from equilibrium to one expo- 
nential in d (uj) (see appendix ICl Theorem 4). 

We want to emphasize that this result about equilibra- 
tion is completely general. That is, we did not assume 
anything special about the interaction (apart from not 
having degenerate energy gaps - which rules out only a 
set of Hamiltonians of measure zero), neither have we 
assumed any special properties of the bath. In particu- 
lar, we did not assume that the bath is characterized by 
some temperature T; in fact we did not assume that the 
bath is in any kind of equilibrium at all! Furthermore, 
we did not make the rather standard assumption that 
the system has a limited spread of energies AE. Finally, 
we also did not make any assumptions about the form 
of the equilibrium state ujs of the subsystem; in particu- 
lar there is no need for this state to have any of the usual 
thermal properties, e.g. to be of Boltzmannian form or 
similar. Indeed, since the bath is completely arbitrary, 
so is the equilibrium state of the system. In other words, 
the equilibration phenomenon that we describe is a gen- 
eral phenomenon, and needs not have any "thermal" as- 
pects at all. 

We also note that the bounds are completely indepen- 
dent from the energy eigenvalues. Indeed, the energy 
eigenvalues did not play any role in establishing the 
bounds - they just gave rise to phases which averaged 
to zero when taking the time average. This is a very im- 
portant feature of our results. In fact, apart from the dis- 
cussion in the second half of section |3] energy does not 
appear anymore. Energy obviously plays a central role 
in establishing the rate of time evolution and its exact 
details, but, as we have shown here, the long-time be- 
havior is characterized by bounds that are independent 
of energy. 

Furthermore, it is not only the energy eigenvalues 
that play no role in our bounds - the form of the en- 
ergy eigenstates is also irrelevant. Indeed, the quantity 

1 / d? 

2 V rfcft'^^) which appears in our theorem is independent 
of the form of the energy eigenstates. 

Finally, in the discussions above we referred to a spe- 
cific partition of our total system into a "subsystem" and 
its corresponding "bath". However, we note that the 
greater upper bound in Theorem 1 depends only on the 
dimension of the subsystem, not on what that particular 
subsystem is, and on the effective dimension of oj which 
is independent of the way we partitioned our total sys- 
tem. Hence when one particular subsystem equilibrates 
then all other subsystems of the same dimension equili- 
brate as well. 
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Equilibration of typical states 

Our theorem described above puts a bound on the 
fluctuation of the state of a subsystem around the time- 
average. The next question is to find the cases in which 
the fluctuations are small so that the subsystem equili- 
brates. As we will now show, almost all quantum states 
have this property. Intuitively the idea is the following. 

Consider a Hilbert subspace Hr of the total Hilbert 
space of large dimension dn- Consider now an arbi- 
trary orthogonal basis in this subspace. First, since the 
dimension of the subspace is large, it is impossible to 
construct all basis vectors from only a few energy eigen- 
states. Thus, generically, each basis vector is a superpo- 
sition of many energy eigenstates. Second a typical state 
1 4') from this subspace will have roughly equal overlap 
with each of the basis vectors. Hence a typical state | 
will have a significant overlap with many energy eigen- 
states, making very large. Consequently, when 
we chose a typical state from a large dimensional sub- 
space Hr we find that the state is such that all the small- 
dimension subsystems equilibrate. 

More precisely, we can prove the following. 

Theorem 2 i) The average effective dimension 

zvhere the average is computed over uniformly random pure 

states 1*) g Hr C H is such that 

(d^'i^))^ > f . (9) 

ii) For a random state j^") e Hr, C H, the probability 
Pr^ {d^^iio) < ^} that is smaller than ^ is ex- 

ponentially small, namely 

Pr* |d<=ff(c^) < ^1 < 2cxp (-cVd^) , (10) 
with a constant c = %^ w IQ-^. 

The proof of (i) is given in appendix lAl and the proof 
of (ii), which makes use of Levy's lemma, is given in 
appendix |Bl 

Point (i) essentially tells us that the average effec- 
tive dimension is larger than half the dimension of the 
Hilbert subspace, so when we draw states from a sub- 
space of large dimension, the effective dimension d'^^{uj) 
of a typical state is large. Point (ii) makes the result even 
sharper, telling us that the probability of having a small 
effective dimension is exponentially small. 

We now apply these results to two specific cases. 

Equilibration of generic states. We first address the 
question of what happens to a generic state of the whole 
system. In other words, what happens to a state chosen 
at random from the total Hilbert space H. The answer to 
this question can be found immediately from the above 
results by taking the subspace Hr, to be the entire space 



H, meaning that in eqs. ^ and ifTOt we replace dR by 
d. Hence, for systems living in high dimensional Hilbert 
spaces, from (|9j and dlOt we expect that the effective di- 
mension (i°*((jj) is of the order of d. Indeed, eq. iflOl l tells 
us that the effective dimension will be smaller than 4 
only in an exponential small number of cases. Conse- 
quently, since d = dsds, the factor djj/d°^{uj) that gov- 
erns the amount of fluctuations will be approximately 
equal to dg/d ^ ds/ds- In the case of systems com- 
posed of large numbers of particles, the dimension of 
the Hilbert spaces of the system and bath grow expo- 
nentially with their number of particles. Hence this ra- 
tio will drop off exponentially with the total number of 
particles whenever the number of particles in the sub- 
system is no more than a constant fraction of the number 
in the bath, and thus each such subsystem equilibrates. 

Equilibration of systems far from equilibrium. In the 

previous subsection we have shown that a typical state 
I '5) is such that any small enough subsystem equili- 
brates. This however is not the end of the story. Indeed, 
one can raise the following legitimate question. Is it the 
case that states in which a subsystem is initially far from 
equilibrium will equilibrate? The point is that states in 
which a subsystem is far from equilibrium are quite rare 
in the Hilbert space so they are by no means "typical". 
Indeed, it can be shown that the overwhelming majority 
of the states in the total Hilbert space are such that ev- 
ery small subsystem is already in a canonical state In 
this section we focus specifically on what happens when 
the initial state of the subsystem is far from equilibrium, 
and show that it too equilibrates. 

We now consider the usual situation in which the 
question of equilibration is discussed. There is a bath, 
consisting of a large number of particles, about which 
we know only some macroscopic parameters (such as 
temperature), into which we place a much smaller sub- 
system. The initial state of the subsystem is arbitrary, 
but uncorrelated with that of the bath. The question is, 
does the subsystem equilibrate? 

In this section we prove that for any initial state of the 
subsystem, and almost all initial states of the bath, the 
subsystem equilibrates. Note that this includes cases in 
which the subsystem is initially far from equilibrium. 

In our formalism, we describe the above situation 
as follows. We consider initial state of the full system 
(i.e. of the subsystem and bath) to be a product state 
|*)sB = |V')s|0)s- The initial state of the system 1-0)5 is 
an arbitrary state in Hs- We model the macroscopic pa- 
rameters of the bath by restricting its initial state 1 0) s to 
lie within a particular subspace ^Hb (of dimension 
dg). Note, however, that although the motivation for 
imposing a restriction on the bath's Hilbert space was 
to model macroscopic parameters, our results apply to 
any restricted Hilbert space Hg. In particular it need 
not have any thermodynamic or macroscopic meaning. 
Furthermore, H§ is only a restriction at the initial time, 
and the state of the bath may evolve outside this space 
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over time. 

Given this setup, we can apply our second theorem, 
eqs. lO and ||lOll with Hr = \ip)s ® / and hence 
Ar = dg. This gives d°^{oj) > dg/A for almost all ini- 
tial states of the bath, and any initial state of the sys- 
tem. Hence given an initial product state, the subsystem 
will equilibrate for any initial state of the system and al- 
most any initial state of the bath (within Ti.^), as long as 
» d|. 

The actual mechanism by which the subsystem equili- 
brates when it is put in contact with the bath is however 
highly nontrivial. Indeed, as we show below, generi- 
cally the state of the bath evolves through many distinct 
states and does not reach equilibrium. Moreover it may 
move out of the initial restricted subspace Hr. At the 
same time, surprisingly, the subsystem that is in contact 
with it does equilibrate. This is even more surprising, 
given the fact that we did not make any assumptions 
as to the form of the interaction between the subsystem 
and the bath (apart from the Hamiltonian not having de- 
generate energy gaps). Thus, in principle, it could be the 
case that the evolution of the subsystem is sensitive to 
the precise state of the bath. 

To show that the bath does not equilibrate generically, 
we will show that d'^^ {lor) the effective dimension of the 
average state of the bath is much larger than d°^{pR {t}), 
the effective dimension of the bath state at any particular 
time, which essentially means that the bath state contin- 
ues to evolve and does not equilibrate to any particular 
state. Indeed, note that because the two systems are in a 
pure entangled state, rank(pB(i)) = rank(p5(t)) < ^5. 
Since the effective dimension of a state is always less 
than its rank, we obtain that 

d'''{pR{t))<ds. (11) 

In contrast, eq. l|A6|l implies that 

rf'=*(c^i3) > rf°^(c^)/ds. (12) 

Hence for a generic state which, as discussed in Theo- 
rem 2 (ii) obeys c?°^(a;) > we have 

d''''{iVR)> — = 7^»ds>d'^''{pR{t)). (13) 



3. INITIAL STATE INDEPENDENCE 

So far, we have addressed one issue associated with 
thermalisation, namely equilibration. In particular we 
have shown that a small subsystem will equilibrate for 
almost all states drawn from a large restricted subspace. 
We now address a second issue - how the equilibrium 
state reached by a subsystem depends on the initial 
state. So far, each initial state obeying the restriction 
could cause the subsystem to equilibrate to a different 
state. To represent this dependency we now explicitly 
denote the equilibrium state of the subsystem by ujg. 



However, when discussing thermalisation, we would 
usually expect the equilibrium state to depend only on 
macroscopic parameters (such as the temperature of the 
bath), and not on the precise microscopic initial state. 
Our next result shows that this intuition caries over to 
our more general setting, by proving that almost all 
states within a large restricted subspace lead to the same 
equilibrium state of a subsystem. 

Theorem 3 i) Almost all initial states chosen from a large 
restricted subspace will yield the same equilibrium state of a 
small subsystem. In particidar, with referring to the av- 
erage over uniformly random pure states |5'(0)) S Hr C 
Hs^Hb, andns - (t^|)^; 




The first inequality gives a tighter hut more complicated 
bound, with 




where Hr is the projector onto Hr. 

ii) for a random state \^) e Hr d H, the probability that 
D{ujg,i}s) > \ + e drops off exponentially with e^dR-. 

Pr* [oiut.ns) > \^|^+^^ < 2cxp {-c'e'dR) , 

(15) 

with a constant c! = Note that setting e = d^^^^ 

yields a small average distance with very high probability 
when dR^ ds. 

The proof is given in appendix |A] [part (i)] and [B] [part 
(ii)]. The quantity 5 reflects the average purity of the 
energy eigenstates on the subsystem (with probability 
weights {Ek\^^\Ek) determined by their relevance to 
H-ji). The quantity 5 lies in the range (^,1), with highly 
entangled energy eigenstates yielding smaller values. 

Bath state independence. We now apply the result 
above to the situation described in the previous sec- 
tion, in which the subsystem and bath are initially in the 
product state \^)sb = \4')s\4>)b in the restricted sub- 
space Hr = IV') ^ H^. Previously, we have shown that 
generic states of the bath cause the subsystem to equili- 
brate. However, the particular equilibrium state of 
the subsystem could in principle depend on the precise 
initial state of the bath 1 0) s . Here we show that this is 
not the case, and that almost all states of the bath in H^ 
lead to the same equilibrium state of the subsystem. 

The proof follows simply by applying Theorem 3 with 
T^R = IV") ^ H^ and hence dR = d^. From the weaker 
bound (not involving 6) we see that almost all initial 
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states of the bath yield approximately the same equilib- 
rium state ojs of the subsystem, as long as ^ ds- 
Note that since we do not use S, this result does not de- 
pend on any properties of the energy eigenstates. 

Subsystem state independence. We now turn to the 
third aspect of thermalisation - subsystem state inde- 
pendence. A cup of coffee left untouched in a closed 
room eventually reaches room temperature regardless 
of its own initial temperature. In other words, a sub- 
system in contact with a bath reaches an equilibrium 
state that only depends on the bath but not on its par- 
ticular initial state. Could we derive this behavior in 
our formalism? The issue of subsystem state indepen- 
dence turns out to be more complicated than the issues 
of equilibration and bath independence and we have not 
been able to completely solve this problem, however we 
present here a number of preliminary results and obser- 
vations that we believe are crucial in understanding this 
question. 

There are many reasons why the issue of subsystem 
state independence is complicated. The first reason is 
that in fact it is not always the case that the equilibrium 
state is independent of the initial state of the subsystem. 
For a drastic example, consider the subsystem to be an 
atomic bomb. Then, depending on whether the bomb 
is functional or not, the room temperature itself will be 
modified, and therefore also the final equilibrium state. 
In other words, although the subsystem appears to be 
small, it may actually have a big impact on the bath, in 
which case the equilibrium state depends on the initial 
state of the subsystem. Here it seems that not only the 
dimension of the Hilbert spaces of the subsystem and 
bath are relevant (as in the previous sections), but the 
value of energy itself. The energy however turns out 
not to be the end of the story, and perhaps not even the 
key aspect, as we show below. 

To start with, one can easily give an explicit example 
of a situation in which the equilibrivun state of a subsys- 
tem depends on its initial state regardless on how big 
the bath is, or on the energy scales involved. In all our 
previous discussions the Hamiltonian was taken to be 
very general. The only condition we imposed on it was 
that it has no degenerate energy gaps. This condition 
is sufficient to ensure that the subsystem interacts with 
the bath ( i.e. H ^ Hs + Hb)- However, this condition is 
not strong enough to ensure that there are no conserved 
quantities of the subsystem. When there are such quan- 
tities, then clearly initial states of the subsystem with 
different values of these quantities cannot equilibrate to 
the same state. An example would be a Hamiltonian of 
the form 

H = ^Enrr^\n){n\s®\m){m\B. (16) 

As long as the energy eigenvalues Enm have no degen- 
erate energy gaps (which could be achieved by choos- 
ing them independently at random from some range). 



this is an interacting Hamiltonian. However, any opera- 
tor A = ^nl"')("'ls the subsystem will commute 
with H , and therefore is conserved. Evolution under 
this Hamiltonian will dephase the state of the subsys- 
tem in the \n) basis, but cannot flip an |n) into an \n'). 

It is tempting to think that we could rule out such 
cases by only considering Hamiltonians that have no 
conserved quantities on the subsystem. However, this 
is still not a sufficient condition, as illustrated by the fol- 
lowing counterexample (in which the system is a single 
spin, and the bath is composed of many spins). 

H = Eal^+Hint + HB, (17) 

where the eigenvalues of -ffint and Hb lie between — 1 
and 1, and E ^ 1. It is easy to choose an Hint such that 
there are no conserved quantities on the subsystem. In 
such cases, one might imagine that in the long time limit 
it would be possible for all initial states of the subsystem 
to equilibrate to the same state. However, it is easy to see 
that is not the case. Consider two initial states |V-'+)s|0) b 
and s\(f>) B where cr^ | ijj± ) = ± | '4'± ) . The difference in 
expected total energy between these two states is close 
to 2E (the contribution of the terms i^int and Hb to the 
expected energy of any state lies between +2 and —2) 
and remains constant in time. If the subsystem were to 
equilibrate to the same state in both cases, this differ- 
ence in energy would have to be accounted for by Hmt 
and Hb- However these terms can generate an energy 
difference of at most 4. Hence such an equilibration pro- 
cess for the subsystem is impossible. 

The above situation is indeed an abstract version of 
the "atomic bomb" example, in the sense that the energy 
of the subsystem completely dominates that of the bath, 
despite the fact that the dimension of the Hilbert space 
of the subsystem is much smaller than that of the bath. 

It is important to note that a similar argument can be 
made for any globally conserved quantity in which the 
total value is dominated by the contribution of the sub- 
system. Hence the issue is not one of energy. In fact 
we believe that the only special role of the energy (more 
precisely, of the energy eigenvalues) is that of fixing the 
time scale of the evolution. Indeed, in all our previous 
arguments the energy eigenvalues drop out of all the 
calculations (as long as there are no degenerate energy 
gaps). Instead, we note that the "atomic bomb" exam- 
ple in its abstract form of the Hamiltonian in eq. (iTl l is in 
fact an approximation of the trivially non-thermalising 
Hamiltonian | [T6l l, as it can be seen quite easily that for 
large E its eigenstates are arbitrarily close to product 
states. 

The above discussion motivates us to consider the 
case in which the energy eigenstates are far from prod- 
uct. In this situation we prove that almost all initial 
states of the subsystem lead to the same time-averaged 
state. Hence when such states do equilibrate (i.e. spend 
almost all times close to their time average) they reach 
the same equilibrium state. 
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We emphasize that as we mentioned in the introduc- 
tion however, the issues of equilibration and subsystem 
independence are in fact totally separated. Indeed, the 
question we address here is that of subsystem indepen- 
dence of the time-averaged state of the subsystem, regard- 
less on the fact whether or not the system has only small 
fluctuations around this state, in which case the time- 
averaged state is an equilibrium state, or not. 

Consider the following situation: let the initial state 
of the system be a product state 1^") = |V')sI'/')b/ where 
now, in contrast to the previous section, 1 0) is fixed but 
IV") G Ti-s is generic. I.e. \'i>) comes from a subspace 
Hb. ^ 'Hs ® \4')b- Again we apply Theorem 3. How- 
ever, we now note that the weaker bound (not involv- 
ing 5) is essentially useless, as dn = ds- Indeed, this is 
to be expected: the weaker bound does not take into ac- 
count any assumptions about the energy eigenstates and 
we have argued above that there are cases in which the 
time-averaged state of the subsystem depends strongly 
on its initial state. 

This leaves us with the stronger bound involving 5. 
Substituting dfj = ds, we see that the time-averaged 
state of the subsystem will be the same for almost all 
initial states as long as 5 ^ 1 . However, this is precisely 
the case we are interested in,in which the relevant en- 
ergy eigenstates are far from product. 

When the energy eigenstates are highly entangled, in 
the sense that (5^1, the equilibrium state of the subsys- 
tem will therefore be approximately equal for almost all 
initial system states. 

4. CONCLUSIONS 

Understanding the basic laws of statistical mechanics 
from first principles is a holy grail of the subject. Here 
we approach this question from a fundamentally differ- 
ent viewpoint than usual, in which the whole system is 
described by a pure quantum state and the probabilistic 
behavior of a subsystem is an objective phenomenon, 
due to quantum entanglement, rather than the result 
of subjective ignorance. We proved a general quantum 
result: With almost full generality all interacting large 
quantum systems evolve in such a way that any small 
subsystem equilibrates, that is, spends almost all time 
extremely close to a particular state. The only conditions 
we require are that the Hamiltonian has no degenerate 
energy gaps (which rules out non-interacting Hamilto- 
nians) and that the state of the whole system contains 
sufficiently many energy eigenstates. 

Virtually all physical situations satisfy these require- 
ments. Firstly, all but a measure zero set of Hamil- 



tonians have non-degenerate energy gaps. Indeed, an 
infinitesimally small random perturbation will lift any 
such degeneracy. Secondly, the vast majority of states in 
the Hilbert space are such that they contain (i.e. have 
significant overlap with) very many energy eigenstates, 
as required in our proof. In particular, this covers the 
physically interesting situation in which a subsystem, 
initially out of equilibrium, interacts with a large bath. 
In this situation, we have proved that for every state 
of the subsystem and almost every state of the bath, 
the subsystem equilibrates. Furthermore, we have also 
proved that the equilibrium state of the subsystem is 
independent of the specific initial state of the bath and 
only depends on its macroscopic parameters. 

We would like to emphasize that the above two re- 
strictions are the only conditions we require. We do not 
require the assumptions usually made in statistical me- 
chanics. For example the interactions could be strong 
and long range and energy need not be an extensive 
quantity. Also the state of the bath could be a super- 
position of a very large range of energy eigenstates and 
hence have no well defined temperature. 

Although our original motivation was to study the 
phenomenon of thermalisation, we have found a much 
larger range of phenomena: we discovered that reach- 
ing equilibrium is an almost universal behavior of large 
quantum systems. Ordinary thermalization is just a par- 
ticular example of this behavior; it occurs in specific sit- 
uations and it has a number of additional characteristics 
that we did not address in the present paper. For ex- 
ample, the equilibrium state of a subsystem that is put 
into contact with a thermal bath, is largely independent 
of the initial state of the subsystem. Also the equilib- 
rium state has a particular form (generally Boltzman- 
nian). These issues are open questions for future work. 

Note added. After the completion of this work we be- 
came aware of an independent work by P. Reimann pub- 
lished very recently, Phys. Rev Lett. 101, 190403 (2008), 
which takes first steps towards solving the thermalisa- 
tion problem along similar lines as in the present paper. 
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APPENDIX A: PROOFS OF EXPECTATION VALUES 

In this appendix, we derive several expectation values 
that are used in the main paper. 

Proof of Theorem 1. We first relate the trace distance 
to a less natural, but more mathematically tractable dis- 
tance measure (the square of the Hilbert-Schmidt dis- 
tance) using a standard bound 



where 



D{pi,p2) = 2 Trs V {Pi - P2Y < -VrfsTrs(pi - pa)^- 

(Al) 

Using the concavity of the square-root function, we 
therefore have 

{Dipsit),u;s)), <^ds (tts (psW-^sf)^- (A2) 

Using the expansions in terms of the energy eigenstates 
given in eqs. ^ and (HJl, 

ps{t) - c^s = E Cfecre-^(^'=-^')* Ttb lE^mi (A3) 

where J2k^i shorthand for a sum over k and I omitting 
all terms in which k = I. Hence, 

{TTs{ps{t)^ujsf)^ 

= J2Y1 Trs(TrB \Ek){Ei\TrB \E^^){En\), 

k^l m^n 



(A4) 



Evaluating this time average, using the fact that the 
Hamiltonian has non-degenerate energy gaps, and that 
the sums only include those terms where k ^ I and 
m ^ n, we find that the only non-zero terms are those 
where k = n and / = rn, giving 

{TTs{ps{t)-iosf), 

= IcfcHc'l' Trs(TrB \Ek){Ei\TTB \Ei){Ek\) 

k^l 

= EI^^-I'l^'l' E {■^b\E,){EiWb){s'y\Ei){E,\sb') 



k^l 



'bb' 



El^^-I'l^'l' E {sb\E,){E,\sb'){s'b'\E,){Ei\s'b) 



k^l 



ss'bb' 



= E \^k?\ci? TrB(Trs \Ek){Ek\i:Ts \Ei){Ei\) 



k^l 



^J2'^TB{TTsi\ck\'\Ekm\)Trs{\cim){E^\)) 

kytl 

= Tib - E I^^I" ((Trs \Ek){Ek\f^ < Ttb luI, 

k 

(A5) 

where ujb = Trs uj. To obtain a further bound, we 
invoke weak subadditivity of the Renyi entropy (see 
e.g. El), 



Tr(.^)> Tr.K) >Tr.K) 



Rank(ps) ds 
and therefore 

{D{ps{t),u^s))t < l\/dsTTB 



(A6) 



<-JdlTr{^^). 



(A7) 



Using the definition of the effective dimension lO, these 
are the bounds as stated in the theorem. □ 

In order to prove Theorems 2(i) and 3(i), we require 
a couple of additional mathematical tools, that are also 
used in |j3|] . The first is the simple and very helpful iden- 
tity 



Tt{AB) = Ti{{A®B)S), 



(A8) 



where S is the SWAP operator of the two systems. This 
equation is easily proved by expanding it in a basis. 
The second result we employ is that (l^-X*! (g) |^')(*|)^ 
is proportional to the projector onto the symmetric sub- 
space. This is actually the simplest instance of the rep- 
resentation theory of the unitary group. It implies the 
following lemma: 
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Lemma For states \^) e Hb. C H 



dR{dR + 1) 



(A9) 



zvhere 1 is the identity operator onH'SQH, and Hrr, = IIj^ ® 
Hr is the projector onto Hr ® Hr. □ 



For conciseness, in the proofs of Theorems 2(i) and 3(i) 
we abbreviate energy eigenstates by their indices, such 
that \Ek) = \k) etc. Similarly, we write \Ek) ® \Ei) as 
|fc^) etc. It is also helpful to note that time-averaging a 
state corresponds to de-phasing it in the energy eigen- 
basis (due to the non-degeneracy of energy levels). We 
denote this dephasing map by %[p\ := \k){k\p\k){k\, 
suchthatc^ = (|V^)(7/.|), = 

With these tools and notation in hand, we now prove 
the results needed for the main paper 

Proof of Theorem 2(1). We first prove a bound on the 
expected purity of lo, 



(TrH2)^ = (Tr((c.®c.)S)>^ 



Tr($(^$[(|*)(*|® |*)(*|)^]S) 

nflfl(i + §)" 



Tr 



dR{dR + 1) 



kl 



E 



dR{dR + 1) 



2{kk\nRR\kk) 

dR{dR + 1) 



< 



^ d 



2{ k\nR\k) ^ 2_ 

R{dR + 1) dR 



(AlO) 



It then follows straightforwardly that 



' ^ >^-^>%, (All) 



Tr(c.)V*- (Tr(c.)2)^ 2 



concluding the proof. 



□ 



Proof of Theorem 3(i). As in the proof of Theorem 1, we 
first relate the trace distance to a less natural, but more 
mathematically tractable distance measure j^, 

(7^(^s,(^s)*))* < ^iv^rf5Tr(^s-(^s>*)'^^ 

< \^ds (Tr(c.s-(c^s)*)')(A12) 



We then prove bounds on the averaged term 

= (TrK))^-Tr5((^5>l) 

= Tr (((tjs ® ^s)^ - (^s)* ® (ws)*) S) 



Trss ^Trssl $ ® $ 



= Trss [Trss 
< Trss Trss 



(|vI/)(^/|®|vI/)(*|)^-IIii®IIii 
dR dR 

ni?j^(l + §) ^RR 



. dR{dR + 1) 



r/2 



s 



d^ 



= ^TrssfTrssf \kl){kl\^\lktkl\ ) § 



^ (fcMfc) T,^(Tr.(|fc)(fc|)TrB(|Oa|); 



4 



{k\^R\l)mR\^ I (Trg \ktk\f + (Tr, 



4 



= 7-E('^l!rl'=>T'-s((T>-B|'n)(t|) 



(A13) 



Note that in the second inequality (9th line) we have 
used the fact that Trs(Tri3 \k){k\ - Trs > (it is 

the trace of a positive operator). 

Inserting these results into eq. (|A12t , specifically, the 
bounds given by the third line from bottom and the last 
line, we obtain the inequalities stated in the theorem. 

□ 



APPENDIX B: USES OF LEVY'S LEMMA 

In this appendix we prove Theorems 2(ii) and 3(ii), in 
which exponential bounds are placed on the proportion 
of states which have values of d {uj) and D{ujg , ils) far 
from the average value, using Levy's Lemma for mea- 
sure concentration on a hyper sphere. 



Levy's Lemma (m Let f : 5^"^ -> R be a real- 
valued function on the {D — l)-dimensional Euclidean sphere 
(which we think of as embedded into D-dimensional Eu- 
clidean space), with Lipschitz constant A = sup^^ l/(^i)~ 
f{x2)\/\xi — X2I2. Then, for a uniformly random point 



Vrx{f{X) > (/)+£} <2cxp 



□ (Bl) 



Note that the pure quantum states G TCr, can be 
thought of as lying on a {2dR — l)-dimensional hyper- 
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sphere, with coordinates given by the real and imag- 
inary components of the state's overlap with an or- 
thonormal basis \h) oiHR: 



X2n-im = M{nm, X2„(*)=lm[(n|vl/)]. (B2) 



In this coordinate system Euclidean and Hilbert space 
norm coincide: - x{^2)\2 = - 1*2) [j- 



Because it involves the simplest application of Levy's 
Lemma, we begin with the proof relating to initial state 
independence. 



Proof of Theorem 3(ii). To prove that almost all ini- 
tial states 1^) yield the same equilibrium state, we ap- 
ply Levy's Lemma directly to the function = 
/(.f(l')) = D{uj'^,ns) on the (2dfl^ - 1) -dimensional 
hypersphere of quantum states. Using an analogue of 
eq. (|C6l l to give 



\D{u^s\ns) ~ D{ul\^ls)\ < ||vl'i)~|vl/2)| (B3) 



we find that the Lipschitz constant of the function satis- 
fies A < 1. Substituting this into Levy's Lemma, as well 
as the average value obtained in Theorem 3(i), we obtain 
the desired result. □ 



Proof of Theorem 2(ii). Unfortunately we carmot prove 
the result in Theorem 2(ii) directly by applying Levy's 
Lemma io d'^^{uj). Instead, we apply it to the function 



/(vI/)^/(x(vI/))=ln Tr $[|vE.)(vl/|] 



(B4) 



where the super-operator $ acts on the subspace TYt Q 
H spanned by energy eigenstates with nonzero projec- 
tion onto Hr (i.e. states |fc) satisfying {k\IlR\k) ^ 0). The 
subspace Ht contains all states which could arise dur- 
ing the evolution of an initial state in Hr, and $ maps 
these states back into H r according to 



~$[p]=Y.\~^){k\p\k){k\ and = _===n^|fc). 

(B5) 

Note that when the Hamiltonian commutes with 11 r, $ 
is identical to the normal $ on Ht- Computing the aver- 



age value of our function we find 
ln(Tr($[|*Xvl/f)))^ 
<ln(Tr($[|vI/)(*|]2))^ 
= InTr ($ ® $ [(|*)(*| ® I^X^Dg,] § 

TlRR{t+S) 



In Tr $ ® $ 



dRidR + 1) 



= In 

\ kl 

< In [ ^(Jfc|M) 



rr.rn f {kl\IlRR{\kl) + \lk)) 



V MdR + 1) 

{kl\URR{\kl) + \lk)) 



(B6) 



kl 



< In 



In 



dR{dR + 1) 



dR{dR + 1) 



dR^df 



To bound the Lipschitz constant of the function 
we actually employ a further function 



g(l') = InTr 



^|n)(n|$[|vl/)(vl/|]|n)(n| 



(B7) 



where \ fi) is an orthonormal basis of Hr. Writing 



tnko = Re[{n\k){k\^)] and t„fei = liji[{n\k) {k\^)] 

(B8) 

it follows that 



.g('J') = InTr 



2 

nkz 



(B9) 

To boimd the Lipschitz constant of g, it is sufficient to 
find an upper bound on its gradient. 

dg 1 

dtnkz ^ y , (y,, ,t^,,, ,y 



.2.2.t,,kzJ2*nk'.'- (BIO) 



Introducing the notation Pn — J2kz ^nkz' noting that 
J2n Pn = 1/ we find 



dg 

d^nkz 



< 



~ (EnPlf 

- 16 

iEnPlf 



(Bll) 
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and hence that the Lipschitz constant of g is at most 

To obtain the Lipschitz constant of /, we note that 
.9(^) > /(^)/ with equality if {\n)} is an eigenbasis of 
Trs 1 I . Now, for any two vectors, we may without 
loss of generality assume that > f{^2), and take 

{|n)} to be the eigenbasis of Tr^ |*2)(*2|- Thus, 



(B12) 

and hence the Lipschitz constant of / is also upper- 
bounded by 4:\/1Ia- 

Applying Levy's Lemma to /(^'), and using the ob- 
servation that Prja; > a} < b and x > y implies 
Pr{y > a} < 6 to substitute the bound on ob- 
tained above, and the further bound 



Theorem 4 Under the above assumption of ergodicity, 



ds 



< exp {-c"eW'*^{uj)) , 



(C3) 



with a constant i 



128^-= 



Proof. We already argued the first equality, so only the 
second inequality is to be proved. 

From Theorem 1 in the main paper, the average of 



In Tr 



> InfTrl 



D{ps{a),uJs) is upper bounded by i 



ds 



, hence 



I |] ) ) (B13) Markov's inequality its median is upper bounded by 



gives 



2e'^ _ f e^y/d^. 



Pivp \ ln(Tr($[|1')(*|]2)) >\n—\< 2 cxp 



727r3 
(B14) 

Negating both sides of the expression inside the braces, 
and then taking their exponent yields 



dn. 



< 2 cxp 



2e^ J ~ V 727r3 

Finally setting e = In 2 yields the desired result. 



(B15) 

□ 



APPENDIX C: FLUCTUATIONS FROM THE 
TIME-AVERAGE 

Here we prove an exponential bound on the propor- 
tion of time the state spends a long way from its equi- 
librium state, using measure concentration results O, 
and a stronger assumption on H. 

Making the assumption that the eigenenergies E of H 
have no rational dependencies (which is much stronger 
than our non-degenerate energy gaps condition), the 
trajectory 



over time fills the torus 



T= ■.= Y,e""'Ck\Ek) : < < 2^| 



(CI) 



(C2) 



uniformly. What this means in fact is that an ergodic the- 
orem holds, equating the time averages of all continuous 
functions of the state with their respective "ensemble" 
averages, i.e. with integration over independent phase 
angles Ofc. The same is true for the indicator functions 
of open sets, as the one in the following theorem. 



ds 



The space on which we want to study the measure 
concentration is the direct product of unit circles. As the 
metric on that space we choose the following weighted 
£^ metric: 



d{a,l3) 



±Y.\c,\^\a,-(3,\. (C4) 



The fc-th factor has diameter \ck\'^, so using El 
Thm. 4.4], we conclude that the concentration function 
i} of this space obeys 

^(r) < e-'-'/*^', where = Y^ \ck\^ = l/d''^{uj). 

k 

(C5) 

"Concentration function" means that for any set A c 
[0; 27r)'' of measure > 1/2, the r-neighbourhood = 
{a s.t. 3/3 e A d{a, [3) < r} has measure at least 1 — t?(r). 

Hence, all that remains is to relate the distance r to the 
variation of the function D{ps{a),uis)- We first obtain 
the general result 

\Dipsia), Lus) - DipsiH), us)] < D{ps{a), psifT)) 

<Dipia),pm 



< \ma))-\i,{P))\^. 



(C6) 



Together with 

|IV'(a))-|V'(/3))| 



< 2Y,\ck\''\ak-Pk\=i7Td{a,/3), 

k 

(C7) 
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this results in 

\D{ps{a),Lus) - D{ps{l),Lus)\ < VI^^d{a,f3). (C8) 
So, to change the value of D{pg,ujs) by more than e 



above the median, one has to go from an element in the 
submedian set a distance of at least r = in our cho- 

47r 

sen metric d. Inserting this into the above concentration 
fimction yields the result. □ 



